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Abstract 

We conjecture that the roots of a degree-n univariate complex 
polynomial are located in a union of n — 1 annuli, each of which is 
centered at a root of the derivative and whose radii depend on higher 
derivatives. We prove the conjecture for the cases of degrees 2 and 3, 
and we report on tests with randomly generated polynomials of higher 
degree. 

We state two other closely related conjectures concerning Newton's 
method. If true, these conjectures imply the existence of a simple, 
rapidly convergent algorithm for finding all roots of a polynomial. 

1 Conjecture concerning annuli 

Let p(z) be a univariate polynomial with coefficients in C. Let zi, . . . , z n be 
its roots. Let (i, . . . , ( n -i be the roots p'. This paper proposes the conjecture 
that zi, . . . , z n lie in a union of n — 1 annuli, one for each of d, . . . , Cn-i- 
The two radii of each annulus are determined from higher derivative values, 
and the inner radius is a constant fraction of the outer radius. The formal 
statement is as follows. 
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Conjecture 1 There exist two universal constants < t\ < 1 < l 2 with 
the following property. Let z\, . . . , z n be the roots of a degree-n complex uni- 
variate polynomial p(z). Let Ci, • • • , Cn-i be the roots of its derivative. Let 
pi, . . . , p n _i be defined by 
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Define annulus 
for j = 1, 



{z : npj < \z-(j\ < i 2 pj} 



,n — l. Then for each i — 1, ... ,n, 
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(1) 

(2) 
(3) 



An example of the conjecture for a particular polynomial is illustrated in 
Fig. 1. 



NOTE ADDED IN REVISION. Recently, M. Giusti, J. Heintz, 
G. Lecerf, L. Pardo, B. Salvy and J.-C. Yakoubsohn have shown (unpub- 
lished communication) that Conjecture 1 is partly true and partly false. In 
particular, they have shown that i\ exists and may be taken to be (y/E—l)/2. 
They have also constructed a family of polynomials showing that i 2 cannot 
exist (i.e., \z — Q\/Pj ~^ °°) as ^ ne degree tends to infinity for this family. 
Their counterexample apparently does not invalidate either Conjecture 2 or 
Conjecture 3 below. A more detailed revision of this paper is forthcoming. 



This conjecture, if true, suggests a simple algorithm based on Newton's 
method for finding all the roots of a complex univariate polynomial. The 
algorithm (as well as two other related conjectures) is described in Section 4. 

Observe that some of the terms in the min in (1) could be infinite, but 
they cannot all be infinite since p^ is a nonzero constant. Observe also that 
the conjecture is invariant under change of variables of the form z' = az + b 
for any nonzero a G C and any b G C. It is also invariant under conjugation 
and rescaling of the polynomial (i.e., replacing p by cp or cp for a nonzero 
cGC). 
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Figure 1: The 'x' marks shown in the figure are the 12 roots of a degree-12 
polynomial plotted in the complex plane. There is a cluster of two very close 
roots and another cluster of three. The '.' markers show the 11 roots of 
the derivative. The j'th derivative root is surround by an annulus with inner 
radius .66pj and outer radius 1.33pj, where pj is given by (1) and the scalars 
.66 and 1.33 are chosen based on experiments in Section 3 as approximations 
to Li and t 2 - 



The formula (1) is motivated by considering the polynomial p(z) = z n — c. 
The roots of this polynomial are the n roots of c. The derivative roots are all 
at 0. For this polynomial, all the terms in (1) are infinite except the k = n 
term, and therefore pj = \c\ l l n for all j. Thus, we could take i\ — t 2 — 1 for 
this restricted class of polynomials. Some additional background is given in 
Section 4. In Section 2 we establish the conjecture for degree 2 and degree 
3 polynomials. In Section 3 we describe a computational experiment to test 
the conjecture for high-degree polynomials. 

Aside from the application to Newton's method, the issue of location of 
roots with respect to the derivative roots is an inherently interesting matter 
that has attracted quite a bit of attention in the literature. A classic result 
in this regard is the Gauss-Lucas theorem, which states that the roots of p' 
lie in the convex hull of the roots of p. The most comprehensive treatment 
of the relationships between the roots of p' and roots of p appears to be the 
monograph of [4]. The present conjecture, however, appears to provide more 
precise information about the location of the roots of p in terms of derivative 
roots than any of the known theorems in [4]. 

2 Polynomials of degree 2 or 3 

It follows immediately from the observations in the previous section that 
the conjecture is true for polynomials of degree 2 since any nondegenerate 
quadratic, after suitable change of variables and rescaling, can be transformed 
to z 2 — 1. Thus, for quadratics, the conjecture is true with l\ — i 2 — 1. The 
conjecture is trivially true in the degenerate case of a quadratic with a double 
root since p± — and Zi — Ci = 0. 

In the case of degree-3 polynomials, first consider the cases of a root with 
multiplicity 2 or 3. For the multiplicity-3 case, the conjecture is true for the 
same reason as in the last paragraph. For the multiplicity-2 case, by rescaling 
and changing variables, we may assume the polynomial is p(z) = z 2 (z — 1). In 
this case, z\ = z 2 = and z% = 1, while Ci — and ( 2 — 2/3. Also, one checks 
that pi = and p 2 = min( v /4/27, (4/27) 1/3 ) « .3849. Since \z 3 - ( 2 \ = 1/3, 

for the conjecture to be true in this case requires an i\ < (l/3)/(y / 4/27) ~ 
0.866 and i 2 = 1. 

The remainder of this section considers the case of a cubic with three 
distinct roots. Again by rescaling, changing variables, and taking complex 
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conjugates, we may assume that the polynomial is monic and that its closest 
pair of roots lie at —1 and 1, and that the third root a lies in Quadrant I. 
Thus, p(z) = (z — l)(z + l)(z — a) where Re a > 0, Ima > 0, and \a — 1| > 2. 
(The latter inequality follows by the assumption that —1 and 1 are closer to 
each other than to a.) One checks that p'(z) = 3z 2 — 2az — 1 hence 

Ci> C2 — • 

Since a 2 + 3 is in the upper half-plane by assumptions made about a, we will 
assume that the branch of square-root appearing in the preceding formula is 
chosen so that \J c? + 3 lies in Quadrant I. We will let Ci be the ' — ' branch 
and ( 2 the '+' branch. 

Let ( u stand for either d or ( 2 - Solving p'(( u ) = yields the identity 
a = (3C - 1)/(2C„) and hence ( u - a = (-£ + l)/(2( u ). 

Let us now consider the three roots {z±, Z2, Z3} = {1, —1, a} in order. For 
z\ — 1, consider the annulus about Ci- By definition, i\ and 12 must be 
universal lower and upper bounds on the quantity 

minl^KCQ/p^CQI^MplCi)! 173 ) 

IC1-1I 

which, after using the facts that p(Ci) = (Ci — 1) (Ci + 1) (Ci — a ) an d Ci — a — 
(-Ci + l)/(2Ci), simplifies to 



( 


(Ci + 1) 2 


1/2 


(d + 1) 2 


l/3\ 




3C? + 1 




2Ci(Ci - 1) 





(4) 



To establish an upper bound on (4), observe that (3C 2 + 1) - 2Ci(Ci - 1) = 
(Ci + l) 2 , hence 

Ki + l| 2 <|3Ci + l| + |2Ci(Ci-l)|. (5) 
Let a* be chosen as the real root of 

(a*) 3 -a*- 1 = (6) 

This a* is close to 1.3247. It follows from (6) that (1 + 1/a*) 1 / 2 and (1 + a*) 1 / 3 
are both equal to a*. 

Turning back to (5), take two subcases depending on the relative sizes 
of the two terms on the right-hand side of (5). Subcase 1 is that \3(f + 
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1| < a* 1 2^i(^i — 1)|. In this case, the right-hand side of (5) is at most 
(1 + a*) ■ |2Ci(Ci — 1)1) and hence the second term of (4) is at most (1 + a*) 1 ^ 3 , 
which is equal to a*. 

Subcase 2 is that \3(f + 1| > a*|2£i(£i — 1)|. In this case, the right-hand 
side of (5) is at most (1 + 1/a*) • |3Ci + 1| and hence the first term of (4) is at 
most (1 + 1/a*) 1 / 2 , which is equal to a*. This establishes the upper bound 
of a* on (4). 

For the lower bound, we use the following cruder argument. Let Qi 
denote the first quadrant. Observe that the function of a given by Ci( a ) — 
(a — y/ a 2 + 3)/3 is analytic in the interior of Qi with a singularity on the 
boundary (at and therefore its real part is harmonic. This means 

that the minimum value of the real part of (1 is attained either for a on 
the boundary of Qj or in the limit for infinitely large a. Observe that £i = 
a(l — Jl — 3/a 2 )/3, and for very large \a\, Jl — 3/a 2 ps 1 — 1.5/a 2 , hence 



£1 ~ l/(2a) which tends to for large \a\. Thus, the minimum real part 
of Ci is attained on the boundary rather than 00. For a on the positive 
imaginary axis (one boundary of Qi), say a = ti, we compute that £i = 
(— V3 — t 2 + ti)/3, which has real part equal to if t > \f?> else real part 
equal to — \/3/3 or greater for t E [0, y/3]. Along the positive real axis (the 
other boundary of Qi), Ci is increasing, as one can check from the derivative, 
hence the minimum value of the real part is again at and is equal to 
Therefore, Re(G + 1) > 1 - v^/3 hence 

|Ci + 1| > 1 - V3/3 w 0.423. (7) 

(This bound would be improved if we also accounted for the constraint that 
\a- 1| > 2.) 

Next we have the following chain of inequalities to analyze the first term 
of (4): 

|3d 2 + l| = |3d 2 + 6Ci + 3-6Ci-6 + 4| 
< |3d 2 + 6Ci + 3| + |6Ci + 6| +4 
= 3|Ci + l| 2 + 6|Ci + l|+4 



= 3|Ci + l| 2 + 6^±^+4 
Isi + L \ 

< [3 + 6/(1 - v/3/3) + 4/(1 - v/3/3) 2 ] - |Ci + 1| 2 , (8) 

where we have used the inequality (7) to obtain the last line. The quantity 
in square brackets is approximately 39.6. 
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The second term can be similarly analyzed: 



|2Ci(Ci 



= |2C 1 2 + 4Ci + 2-6Ci-6 + 4| 

< 2|Ci + l| 2 + 6|Ci + l|+4 

< [2 + 6/(l-v / 3/3)+4/(l-v / 3/3) 2 ]-|C + l| 



(9) 



Combining (8) and (9) establishes a rather poor lower bound of ^1/39.6 on 



(4). 



For Z2 — — 1, we must obtain lower and upper bounds on 



m\n(\2 V {Ci)/p'\Ci)\ 1/2 Ap(Ci)\ l/ ") 
lCi + l| 



which simplifies to 
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l/3\ 
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2&(& + l) 





(10) 



To establish a upper bound on (10), use an argument analogous to the pre- 
ceding analysis of z\. (3( 2 + 1) — 2£i(Ci + 1) = (Ci — l) 2 , hence 

ICi-i| 2 < |3Ci 2 + l| + |2Ci(Ci + l)|. 

Depending on the relative sizes of the terms on the right-hand side, either 
ICi - 1| 2 < (1 + 1/a*) • |3Ci 2 + 1|, implying that the first term of (10) is at 
most a*, or else |(i — 1| 2 < (1 + a*) ■ |2^i (Ci + l)l> implying that the second 
term of (4) is at most a*. Here, a* was defined by (6). 

To establish a lower bound, again similar arguments are used. We first 
claim that ReCi < 0; this follows again by considering the extremal cases for 
ReCi as a function of a as above. This implies that \d — 1| > 1. 

Then we have the following chain of inequalities to analyze the first term 
of (10): 



|3C 2 + 1| = |3C 2 -6Ci + 3 + 6Ci-6 + 4| 
|3C 2 -6Ci + 3| + |6Ci-6| + 4 



< 



3|Ci 
3|Ci 



l| 2 + 6|Ci-l|+4 



II 2 + 6 



Ki-l| 



10 — 11 

< (3 + 6 + 4)|Ci-l| 2 , 



+ 4 



7 



where we have used the inequality |Ci — 1| > 1 to obtain the last line. 
The second term can be similarly analyzed: 

|2Ci(Ci + l)| = |2d 2 -4Ci + 2 + 6Ci-6 + 4| 

< 2|Ci-l| 2 + 6|Ci-l|+4 

< (2 + 6 + 4)|C- 1| 2 . 

The last root to analyze is z% = a; for this root we will consider the 
annulus about derivative root (2 instead of Ci- The quantity to analyze is 

min(|2p(C 2 )//(C2)| 1/2 ,b(C 2 )| 1/3 ) 
IC2 - a\ 

which, after simplification, is equal to 



mm 
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4CI 


1/2 


4CI 


l/3\ 




3d + 1 


1 


CI - 1 





(11) 



An upper bound on (11) is obtained by observing that 4£f = 3£f + 1 + Q — 1, 
hence 

|4C 2 2 |<|3Cl + l| + lC 2 2 -l|- 
Then, using the same logic as in the previous two cases, we conclude that 
(11) has a* as an upper bound. 

For the lower bound, observe that ^/ a 3 + 3 lies in Quadrant I (denoted 
Qi) provided that a E Qi- If w\, u> 2 are any two complex numbers both lying 
in Qj, then \w\ + w 2 \ > Therefore, since C2 = (o + V^ 3 + 3)/3, we 

conclude that \( 2 \ > |a|/3. The assumptions \a — 1| > 2 and a E Qi together 
imply that \a\ > \/3, and therefore |^| > n/3/3, hence |£f | > 1/3. 

The inequality derived in the last paragraph yields a lower bound on both 
terms of (11). For the first term, |3C 2 2 + 1| < 3|C 2 2 | + 1 < 6|C 2 2 |. Therefore, the 
first term of (11) is at least (2/3) 1 / 2 . For the second term, a similar use of 
the previous paragraph shows \Q — 1| < 4|£f |, and hence the second term is 
at least 1. 

This concludes the analysis of the n = 3 case. We have shown an upper 



bound of a* for t 2 and a lower bound of yl/39.6 for l±. We have written 
a Matlab program that computes i\ and t 2 for each cubic polynomial with 
roots at —1,1, a, where a ranges over a fairly dense grid lying in the set 
{a E Qi : |a — lj > 2}. We found that i\ appears to be approximately 0.82 
while l 2 appears to be exactly a*, and in particular, the upper bound of a* 
on (11) appears to be tight. 
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3 Computational experiment with higher de- 
gree polynomials 

In this section we describe our Matlab computational experiment with higher 
degree polynomials. 

We experiment with three degrees: n = 10, 20, 40. For each n, we generate 
3000 random polynomials. Each polynomial is chosen by selecting its n roots 
uniformly at random on the unit circle. The rationale for this choice (as 
opposed to a distribution over a 2-dimensional domain) is to greatly increase 
the likelihood of nearby or clustered roots, which is a more difficult case for 
root-finding. 

The roots of the derivative polynomial are then computed, as are the pa- 
rameters pi, . . . , p n -i given by (1). For each polynomial root Zi, i = 1, . . . , n, 
the program seeks the j in {1, . . . , n — 1} such that \z^ — Cj\l Pj is closest to 1. 
Call this quotient i Zi . The program then tabulates the minimum and maxi- 
mum i z . encountered among all 3000 trials; these are taken to be estimates 
for i\ and i 2 . 

The results £1X6 clS follows. For n = 10, L\ = 0.67 and 12 = 1.32; for 
n = 20, Li = 0.66 and 12 = 1.33; for n = 40, L\ = 0.66 and 12 = 1.33. Thus, 
there seems to be little appreciable change in the experimental values of Li 
or l 2 as the degree increases. 

We mention two subtleties concerning the implementation of this compu- 
tational experiment. As mentioned above, the program selects the roots of 
the polynomial at random on the unit circle and then computes the deriva- 
tive roots. The naive method to compute derivative roots, namely, form the 
standard monomial basis for the polynomial, differentiate it term by term, 
and then use the Matlab roots function on the derivative, is unstable for 
polynomials with clusters of roots. This naive implementation gave incorrect 
experimental results. We found that a better method for finding derivative 
roots is to compute the eigenvalues of the (n — 1) x (n — 1) matrix 

M = diag(z 2 , ...,z n )- ev T /n 

where e is the vector of all l's and v is the vector whose ith entry is z i+ i — Z\. 
A brief explanation of why these eigenvalues are derivative roots is as follows. 
If Mx = £x, then ZiXi — d = (xi for each % = 2, . . . , n, where d = v T x/n. 
Solving yields X{ = dj {zi—Q. Substituting this formula for X{ into d = v T x/n 
and simplifying yields — zi) + ■ ■■ + !/(( — z n ) = 0, which is the same 
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as p'(C) = 0. 

It is easy to check that this method for derivative roots works much 
better than the naive method for contrived examples of polynomials with 
root clusters, e.g., the polynomial (z — 2) 10 — .01 10 = 0, which has 10 roots 
lying on a circle of radius 0.01 about the point 2. We have not, however, 
attempted a formal proof of stability of this method. 

A second stability subtlety is the computation of pj. The naive method, 
namely, forming p in the standard monomial basis and differentiating term 
by term to obtain all the derivatives in (1), is unstable. Our implementation 
uses the following method, which appears to be stabler. Form the standard 
monomial representation of the polynomial q(z) = p(z + Q) by multiplying 
together its degree-1 factors (z — Z\ + Q), . . . (z — z n + Q) (or equivalently, 
by applying the Matlab poly function to the n shifted roots z\ — Q, . . . , 
z n — Q). Then derivatives of the form p^ k \Cj) are directly obtained from the 
coefficients of q. 



4 Application to Newton's method 

Recall that Newton's method for finding a root of a complex function is 
given by the iteration 

x k+1 = x k -(P(x k )/(P\x k ). 

Newton's method is known to converge quadratically to a nondegenerate root 
x* if the starting point x° is sufficiently close to x*. The basin of attraction 
for x* is the set of starting points x° such that Newton's method will con- 
verge to x* for that starting point. Note that although Newton's method is 
asymptotically quadratic, a point in the basin of attraction could lead to a 
sequence of iterates that meanders far away from x* for an arbitrary number 
of iterations before asymptotic quadratic convergence takes hold. 

This difficulty leads us to define the basin of fast convergence for x* to 
be the set of x° such that the sequence of iterates generated by Newton's 
method starting from x° converges quadratically immediately (rather than 
asymptotically) according to the following inequality: 

\x k -x*\ < Ul/2) \x°-x*\ 
for all k > 0. This definition is similar to one from Blum et al. [1]. 
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For a degree-n polynomial p such that the roots of p' are Ci, . . . >Cn-i> 
let us define the DR-circles to be the set of circles about Ci>--->Cn-i °f 
radius pi, ■ ■ ■ , p n -i respectively. (Here, DR stands for "derivative root.") By 
Conjecture 1, the roots of p apparently all lie close to the union of its DR- 
circles. This suggests that many points on DR-circles lie in the basins of fast 
convergence of the roots of p. 

There are several different possible conjectures that could be made about 
DR-circles and basins of convergence. The first conjecture is that the basin 
of fast convergence of each root of p contains a subsegment of at least one 
DR circle. Consideration of the polynomial z n — c indicates that the length 
of this segment could be as small as 0(1/ n) radians. This is our conjecture: 

Conjecture 2 There is a universal constant r\\ > with the following prop- 
erty. Let p be a degree-n univariate complex polynomial whose roots are 
z±, . . . , z n . Let Ci, . . . , C n -i be the DR-circles of p, that is, circles centered 
about the roots Ci, • • • , Cn-i of p' such that the radius of Cj is pj as defined 
by (1). Let zi be any root of p of multiplicity 1. Then there exists a segment 
of a DR-circle of length rji/n radians lying in the basin of fast convergence 
of Zi. 

To test this conjecture, we used the same set-up (3000 polynomials of 
degrees 10, 20 and 40) as in Section 3. We discretized each DR-circle with lOn 
evenly spaced points, and for each root of p and each DR-circle we counted 
the number of such DR-circle points in the root's basin of fast convergence. 
If the conjecture were true, this number would always be greater than a 
positive constant 10r/i/(27r) for at least one DR-circle. In fact the minimum 
number for n — 10 was 7, for n = 20 was 10, and for n = 40 was 13. This 
gives some evidence in favor of the conjecture. 

If this conjecture were true, it would imply a very simple algorithm based 
solely on Newton's method for finding all the roots of a degree-n polynomial. 
First, find the unique root of the linear polynomial p^ l \ From this root 
of p( n_1 ), find the two roots of p( n_2 ) by starting Newton's method from a 
sufficient number of sample points on the DR-circle of p( n ~ 2 \ Once the two 
roots of p( n ™ 2 ) are found, construct the DR-circles of p( n ~ z \ sample them 
with points, and carry out Newton's method to find roots of p^ n ~ z \ etc., 
until finally we find the roots of p from those of p' . The conjecture suggests 
that the number of sample points per DR-circle ought to be 0(n). 

The complexity of this algorithm may be estimated as follows. Suppose 
the roots lie in a disk of radius R and root accuracy of e is desired. Starting 
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from the basin of fast convergence, Newton's method requires 0(log \og(R/ e)) 
iterations to achieve the desired accuracy. See Renegar [5] for a more careful 
explanation of the factor loglog(-R/e) as well as a matching lower bound. 
For finding the roots of p( n ~ k \ we require Newton's method to be started on 
k — 1 circles, with 0(k) points per circle. Each iteration of Newton's method 
requires O(k) arithmetic operations. Thus, the number of operations for the 
roots of p (ra ~ fe ) is 0(k 3 loglog(i?/e)). This is summed from k = l,...,n, 
yielding a bound of 0(n 4 log log(i?/e)) operations. 

The previous literature has many algorithms for finding all roots of a 
univariate polynomial; see e.g., the survey of Pan [3]. Our complexity bound 
is worse than published bounds for rootfinding in terms of its dependence on 
n, although it is much simpler than most algorithms. 

Another rootfinding algorithm that uses only Newton's method is due to 
Hubbard et al. [2] and is even simpler than ours in that it uses a fixed set of 
Newton starting points that depends only on the degree n. The drawback 
of the algorithm of Hubbard et al. is that Newton's method is in general not 
necessarily initiated in the basin of fast convergence, so complexity estimates 
are far from optimal. 

Renegar also has an algorithm [5] for all roots of a univariate polynomial 
based primarily on Newton's method. Renegar's algorithm always initiates 
Newton's method in the basin of fast convergence and hence also has a run- 
ning time proportional to log log(i?/e) but has a better dependence on n than 
ours. Renegar's algorithm, which to some extent motivated the present work, 
is based on the following key idea. A root z of a polynomial p has a large 
basin of fast convergence unless z is part of a cluster of closely spaced roots. 
Suppose, for example, that three roots of p are clustered, and no other root 
is nearby. In this case, all three will have small basins of fast convergence. 
On the other hand, it is guaranteed in this case that there is a root of p" 
close to the three clustered roots, and this root of p" will have a large basin 
of fast convergence. More generally, an isolated cluster of k roots must be 
near a root of p( fc_1 ) that has a large basin of fast convergence. Thus, Rene- 
gar's algorithm consists of zooming in on root clusters (possibly recursively, 
since clusters can be nested inside other clusters) by using Newton's method 
for roots of derivatives. Once a point near the cluster is found, subdivision 
is used to find the basin of fast convergence for each individual root in the 
cluster. The drawback of Renegar's method is that, in addition to Newton's 
method, it involves some other operations such as computation of approxi- 
mate winding numbers that might be difficult to implement in practice. 
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A proof of Conjecture 2 might lead to further insight that would re- 
duce the n 4 factor in the complexity bound. For example, suppose it were 
possible to predict which DR-circle would have at least one constant-sized 
segment in the union of basins of fast convergence. In this case, we could 
modify the above procedure by tracking only one DR-circle per derivative 
and sampling that circle with a constant number of points. This yields an 
0(n 2 log log(i?/e)) algorithm for finding a single root of p. Then this root 
could be used to deflate the polynomial, and the process could be repeated, 
yielding an 0(n 3 log log(i?/e)) algorithm to find all roots. 

We check this latter possibility by computing, for each polynomial in our 
test set, what is the minimum number of DR-circles that have at least 1/10 
of their sample points (i.e., a total segment length of 27r/10 radians) in a 
basin of fast convergence. This number appears to grow linearly: for n = 10, 
there were always at least 6 such DR-circles, for n = 20 there were always at 
least 13, and for n = 40 there were at always at least 25. Let us state this 
as another conjecture. 

Conjecture 3 There exist two universal constants r\2 > 0, 773 > as follows. 
For any degree-n polynomial p, at least 7/ 2 n of its DR-circles contain segments 
of length i] 3 radians that lie in the union of the basins of fast convergence of 
roots of p. 

A final possible conjecture concerns the total length of basins of conver- 
gence. The total length of DR-circles in radians is 2ir(n — 1), and p has n 
roots, so one might conjecture that the basin of fast convergence of a par- 
ticular root meets DR-circles in at least a constant number of total radians. 
Our computational experiment, however, did not support this conjecture — 
in fact, our experiments suggests that it is more likely that the minimum 
number of total radians in the basin of convergence for a particular root is 
0(l/n), i.e., no better than what Conjecture 2 implies for a single DR-circle. 

Again, there is one subtle numerical stability issue concerning the tests 
in this section. If the Newton update term p(z)/p'(z) is computed by writ- 
ing down p and p' in standard monomial form and substituting the current 
iterate, then an unstable procedure results and the computational test yields 
false results. Instead, the computational test uses the mathematically equiv- 
alent formula 
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This formula is applicable only in the case that the roots of p are already 
known, which obviously would not occur in practical application of Newton's 
method. 

5 Concluding remarks 

This paper raises three conjectures concerning the location of roots of a 
polynomial. Clearly, the main topic for future work would be to prove one of 
them. Another interesting question for practical work concerns a numerically 
stable implementation of the root-finding procedure outlined in the previous 
section. We imagine that for numerical stability, Newton's method should 
always be applied to shifted polynomials (i.e., polynomials of the form q(z) = 
p(z — s), where p is the original polynomial), where the shift is close to the 
sought-after root. 

A final interesting question is whether the conjecture generalizes to arbi- 
trary entire functions whose roots have finite multiplicities. In this case, k 
appearing in (1) ought to run from 2 to oo. 
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